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16 Abstract 

A simple theoretical model which can be used as an aid in the design of the baffle aperture 
region of a hollow cathode equipped ion thruster is developed. An analysis of the ion and elec- 
tron currents in both the main and cathode discharge chambers is also presented. From this 
analysis a model of current flow through the aperture, which is required as an input to the design 
model, is developed. This model is verified experimentally. The dominant force driving electrons 
through the aperture is shown to be the force due to the electrical potential gradient.. The 
diffusion process is modeled according to the Bohm diffusion theory. A number of simplifications 
are made to limit the amount of detailed plasma information required as input to the model to 
facilitate the use of the model in thruster design. This simplified model is shown to give re- 
markably consistant results with experimental results obtained with a given thruster geometry over 
substantial changes in operating conditions. The model does, however, appear to be uncertain to 
about a factor of two for different thruster cathode region geometries. The design usefulness is 
limited by this factor of two uncertainty and by the accuracy to which the plasma parameters 
required as inputs to the model can be specified. 
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I . INTRODUCTION 


In the past the design and optimization of ion thrusters has been 
accomplished primarily by the cut and try method .^ -6 That is, a 
systematic variation of the parameters that influence thruster opera- 
tion was performed until an acceptable configuration was obtained. 

This procedure is both costly and time consuming. Consequently, there 
is a need for the development of analytical models which relate various 
aspects of thruster operation to, a priori, specifiable thruster 
parameters. The present investigation is an attempt at the development 
of a model which deals with the creation of high energy electrons within 
the discharge chamber of a hollow cathode equipped ion thruster. 

After presenting a brief qualitative description of the operation 
of electron bombardment ion thrusters, the reasons for the use of hollow 
cathodes as well as the special considerations required by their use 
will be given. Next, the development of the theoretical model and the 
assumptions used in its derivation is presented. Finally, a comparison 
of the model with experimental data will be given for a wide variety of 
thruster operating conditions and two different discharge chamber 
geometries. 

The plasma properties and operating parameters given are repre- 
sentative of mercury which was the propellant used in this investigation. 
Units will be SI unless otherwise stated. 



II. THRUSTER OPERATION 


Electron Bombardment Thrusters 

An electron bombardment ion thruster typically consists of a 
cylindrical discharge chamber bounded on the sides by an anode, on the 
upstream end by the body, and on the downstream end by the accelerator 
system as suggested by Fig. 1. An electron emitter or cathode is lo- 
cated inside the thruster and a plasma discharge is sustained between 
the cathode and an electron collector or anode which is typically biased 
30 to 40 volts positive of the cathode. The accelerator system consists 
of a set of multi-aperture electrode plates called grids. The inner 
electrode is called the screen grid and is usually maintained at a high 
positive potential, as is the thruster body. The outer electrode is 
biased strongly negative of the screen grid and is called the accelerator 
grid. A strong electric field exists between the grids which acceler- 
ates the ions produced in the plasma to provide the thrust. 

During operation neutral propellant gas is injected into the 
thruster at the upstream end of the discharge chamber. The cathode 
which is heated to incandescence emits a stream of electrons. Since 
the cathode is maintained 30 to 40 volts negative of the anode these 
electrons are accelerated into the main discharge region. Electrons 
which have undergone this acceleration are called primary electrons. 

The energy of these primary electrons is sufficient to assure' that the 
associated ionization cross section is near a maximum for the given pro- 
pellant. A second group of electrons originates from the inelastic 
interaction of the primaries with the ions and neutral propellant 
atoms. These interactions reduce the energy of the primary electrons 
and release low energy secondary electrons to form an electron 
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Fig. 1 Electron Bombardment Ion Thruster 
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population with a nearly Maxwellian energy distribution. The Max- 
wellian group typically contains 90-95 percent of the plasma electrons 
yet only accounts for roughly half of the ionizations. ^ It is possible 
for both electron populations to exist simultaneously in the discharge 

chamber due to the low interaction rate between the primary and 

8 

Maxwellian electrons. 

Since the mean free path for electron-neutral ionization collisions 

is on the order of meters while typical discharge chamber dimensions 

are on the order of several centimeters a magnetic field is employed to 

prevent the primary electrons from having direct access to the anode. 

The strength of the magnetic field is determined by the need to contain 

the primary electrons while at the same time allowing the Maxwellian 

electrons to diffuse to the anode at a rate sufficient to sustain the 

discharge. The configuration and strength of this magnetic field has a 

significant effect on the operation of this type of thruster and con- 

1_4 9-12 

sequently has been the subject of numerous studies. * 

The plasma produced within the discharge chamber will typically 

13 

assume a potential slightly positive of the anode. Consequently a 
potential sheath will exist at all plasma boundaries. The magnitude of 
this sheath at cathode potential surfaces depends on the magnitude of 
the discharge voltage applied between the cathode and the anode, and 
is typically 30 to 40 volts. This sheath potential is sufficient to 
reflect all but the most energetic electrons in the tail of the Max- 
wellian distribution. ^ At the anode, however, there is a significant 
probability that most electrons will be collected since the plasma 
potential is within a few volts of the anode potential . Consequently, 
the vast majority of electrons can only leave the discharge chamber at 
anode potential surfaces. 
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The ions produced within the plasma move with nearly equal prob- 
ability in all directions. Those that fall on solid thruster surfaces 
recombine with electrons and are vaporized back into the discharge 
chamber as neutral atoms. ^ Those that reach the downstream end of the 
thruster are extracted out through the accelerator system to form the 
beam current. 

Electron Sources 

Steady state thruster operation requires that high energy elec- 
trons be continuously supplied to the main discharge chamber to sustain 
the discharge while low energy electrons are continuously removed by 

the anode. Early ion thruster designs used refractory metal filaments 

1 4 

as the source of primary electrons. These filaments, typically made 
of tantalum or tungsten, emit electrons when heated to sufficiently 
high temperatures. The emitted electrons are then accelerated across 
the cathode potential sheath which surrounds the filament into the main 

discharge chamber to become primary electrons. Since the potential 

1 5 

sheath is electrostatic in nature, the acceleration of the electrons 
from the cathode filament is also an electrostatic process, and the 
primary electrons enter the main discharge region with a nearly mono- 
energetic energy distribution. This type of cathode however proved 
unsuitable for space applications due to its high specific heating 
power characteristic (670 W/A) and short lifetime (^100 hours). ^ In 
order to reduce the specific heating power, various cathode designs 
utilizing oxide coatings have been investigated. The oxide coating 
serves to reduce the work function of the emitter allowing it to operate 
at lower temperatures and hence lower specific heating powers. While 
the specific heating power was generally low enough for most oxide 
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cathodes they had the disadvantage that exposure to air caused chemical 

deactivation of the oxide which prevented them from being preflight 

1 4 

tested in a flight qualified system. 

The hollow cathode has replaced both refractory metal and oxide 
cathodes as the source of electrons in ion thrusters for space applica- 
tions for the following reasons. First, they can be operated at specific 

1 4 

heating powers of less than 40 W/A. Second, they exhibit very long 
lifetimes, ' ® (lifetimes in excess of 10,000 hours have been 
achieved).^ This is necessary since the low-thrust nature of electron 
bombardment thrusters requires long mission times to achieve the desired 
spacecraft velocity increment. Finally, a hollow cathode can be restarted 
a number of times even after exposure to air, allowing it to be pre- 
flight tested.^’ 

Hollow Cathode Equipped Thruster Operation 

A typical hollow cathode equipped ion thruster is shown schemati- 
cally in Fig. 2. The use of a hollow cathode generally requires the 
addition of a small loop anode called a keeper located downstream of 
the cathode. This keeper electrode draws sufficient electron current 
to maintain the cathode discharge in the presence of low emission con- 
ditions and plasma fluctuations. When heated to sufficiently high 
temperatures the hollow cathode will emit a stream of electrons. These 
electrons flow out of a small orifice at the downstream end of the 
cathode and into a small subchamber called the cathode discharge region. 
This region is separated from the main discharge chamber through the 
use of a baffle plate. The plasma potential in the cathode discharge 
region will typically assume a potential slightly positive of the 
keeper electrode. Therefore, biasing the keeper positive of the cathode 
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Fig. 2 Hollow Cathode Equipped Electron Bombardment Ion Thruster 


8 


will create a favorable potential sheath near the cathode orifice which 
facilitates the extraction of electrons from the cathode. The plasma 
potential of the cathode discharge region is generally about 10 volts 
positive of the cathode itself. Since the boundaries of this region 
are maintained at cathode potential a potential sheath will exist at 
each of these surfaces. The magnitude of this sheath is sufficient to 
reflect nearly every electron incident upon it. The cathode chamber 
connects to the main discharge chamber via an annular opening called 
the baffle aperture. Because the sheath potentials in the cathode 
discharge region prevent the electrons from reaching the cathode cham- 
ber walls, electrons may only leave this region either by being collect- 
ed by the keeper electrode or by diffusing through the baffle aperture. 

Cathode Discharge Region Current Balance 

Fig. 3 indicates the currents of electrons and ions into and out of 
the cathode discharge region of a normally operating thruster, The ar- 
rows for both ions and electrons point in the direction of actual parti- 
cle motion. Consequently, for the electrons the arrows point in the 
direction opposite to that of conventional current flow while for the 
ions they point in the conventional current flow direction. Linder steady 
state operating conditions the net rate at which electrons enter the 
cathode discharge region must be equal to the net rate at which ions 
enter, that is, 



Electrons Electrons Ions Ions 

Entering Leaving Entering Leaving 


As stated above the electron current to the walls (I ew ) ma y be neg- 
lected. In addition the ion current which enters through the baffle 
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El ectrons 


I = Cathode Electron Emission Current 
e E 

I = Baffle Aperture Electron Current 

I = Electron Current to Cathode Potential Surfaces 
e w 

I = Keeper Electron Current 
e k 

I = Net Cathode Emission Current, I,. = I - I 
E E e E e k 

I = Electron Production Rate 
6 P 

Ions 

I. = Baffle Aperture Ion Current 

v I. = Ion Production Rate 

n p 

I. = Ion Current to Cathode Region Walls 
1 w 


Fig. 3 Electron and Ion Currents in Cathode Discharge Chamber 
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aperture is negligible primarily as a result of the small area of the 
* 

aperture. Further, the production rate of ions (I,* ) and of electrons 

P 

(I 6 p) in this region through ionizations must be equal, therefore Eq. (1) 
may be written in the form. 


I = I - I + I. 
e a e E e k 7 w 


( 2 ) 


Defining the net cathode emission current (1^) as the electron current 
which is emitted by the cathode and not subsequently collected by the 
keeper, i.e. as, 


h = l e " l e 
11 e E e k 


Eq. (2) can be written. 


I = I_ + I. . 
e a E ’w 


(3) 


(4) 


In addition, the calculations presented in Appendix A based on measured 

cathode region plasma densities indicate that the ion current to the 

cathode chamber walls (I,* ) is small compared to the net emission cur- 

w 

rent (1^). Thus Eq. (4) may be written in the form, 

I = Ir .' (5) 


That is, the magnitude of the electron current through the baffle 
aperture into the main discharge region is approximately equal to the 
net cathode emission current. 

The cathode discharge region electrons are driven through the 
baffle aperture under the influence of a density gradient force and a 


Typical ion currents through the aperture are on the order of 0.05 
amps, and are less generally than 2% of the electron current 
through the aperature. See Appendix A. 
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potential gradient force. The density gradient arises because the 
hollow cathode creates a cathode region plasma density which is gen- 
erally greater than the density of the main discharge region. The 
potential gradient occurs as a result of the anode being biased posi- 
tive of the cathode. A magnetic field oriented across the aperture is 
required to control the impedence of the baffle aperture region, 3,4, 6,1 7-21 
This creates the following self consistant situation. The anode, 
being maintained at a potential positive of the cathode, draws the elec- 
trons emitted by the cathode through the baffle aperture into the main 
discharge region. The magnetic field oriented across the aperture cre- 
ates an impedence to this flow of electrons, and thereby establishes 
the potential and density gradients in the vicinity of the baffle aper- 
ture. These potential and density gradients then serve to drive the 
electron current across the magnetic field and into the main discharge 
region. Hollow cathode equipped ion thrusters typically operate with a 
potential increase of about 30 volts through the baffle aperture. Elec- 
trons diffusing through the aperture region thus acquire an energy of 
approximately 30 electron volts and subsequently are referred to as 
primary electrons. In this case, however, the accelerating potential 
sheath is electromagnetic in nature,^ 7 not electrostatic as it was for 
the refractory metal and oxide cathodes. These primary electrons now 
have sufficient energies to ionize the neutral propellant atoms. 

Main Discharge Region Current Balance 

Electrons may enter the main discharge region either by coming 
through the baffle aperture or through the liberation of electrons in 
the ionization process. In both cases the electrons, as stated pre- 
viously, are constrained by the plasma boundary potential sheaths to 
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leave the main discharge region only at the anode. The ions on the 
other hand only enter the main discharge region through ionizations. 
They may, however, leave this region by going through the baffle 
aperture, recombining on a cathode potential surface, recombining on an 
anode potential surface or being extracted into the ion beam by the ac- 
celerator system. Again for a steady state operating condition the net 
influx of electrons must be equal to the net influx of ions. With the 
help of Fig. 4, where the arrow convention is the same as that used in 
Fig. 3, the following equation can be written. 


I + I - I -I =1. - I. 

e a e P e w e A ’P ^A_ 


I. - I. - I. 
’w ’a ’b 


( 6 ) 


Electrons Electrons Ions 
Entering Leaving Entering 


Ions 

Leaving 


Now since the electron current to cathode potential surfaces (Ig^) anc * 

the ion current through the baffle aperture (Ij ) are negligible and 

since the production rate of electrons (I e ) and of ions (I.- 1 through 

r 'P 

ionizations are approximately equal (neglecting double ion production) 
Eq. (6) can be written in the form. 


I 


I 


- I. - I. -I. 


( 7 ) 


a A 'A 'w 'b 

Because both electrons and ions are collected at the anode the dis- 
charge current, Ip, (i.e., the current through the anode power supply) 
is simply the difference between the electron and ion currents to the 
anode. 




(8) 


Combining Eqs. (7) and (8) yields, 


I = I n - I. -I. 
e D i,, i k 
a w b 


( 9 ) 
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El ectrons 

I = Baffle Aperture Electron Current 

I e M Electron Production Rate 

I » Electron Current to the Anode 
A 

I g = Electron Current to Cathode Potential Surfaces 
w 

Ip = Electron Current Through Discharge Supply 
I s = Electron Current Through H V Screen Supply 
= Net Cathode Emission Current 


Ions 

I. = Baffle Aperture Ion Current 

I = Ion Production Rate 
P 

I = Ion Current to the Anode 

’a 

I. * Ion Current to Cathode Potential surfaces 
’w 

I. * Ion Beam Current 
'b 


Fig. 4 Electron and Ion Currents in Main Discharge Chamber 
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Assuming the ions move with nearly equal probability in all directions 

the ion current to the cathode potential surfaces may be determined in 

terms of the beam current (I n -, ) by the relation, 

1 b 



where A $ is the effective open area of the screen grid and A w is the 

★ 

effective surface area of the ion production region through which ions 
can flow to reach a cathode potential surface. Combining Eqs. (9) and 
(10) yields. 


By measuring the discharge current and beam current one should be able 
to approximate the current through the baffle aperture using this 
equation. Equation (11) thus provides an alternate method for approxi- 
mating the baffle aperture current when it is not possible to measure 
the net cathode emission current directly. 

Present Investigation 

The present investigation is an attempt at the development of a 
theoretical model which relates, a priori, specifiable magnetic field, 
geometric and plasma parameters (density, potential, etc.) to the 
aperture impedence. It was not the intent of this investigation to 
determine the exact mechanism of the diffusion process of electrons 
through the aperture. Rather the simplest relation between the 



* 


The ion production region is defined as the region of the main 

13 

discharge chamber inside the virtual anode surface. 
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electron current through the aperture and the forces acting on-it, which 
agreed with experimental results, was sought. As stated previously, it 
is believed that the current is driven through the aperture across the 
magnetic field under the influence of both a plasma density gradient 
and a potential gradient. Thus any model describing the flow of elec- 
tons in this region would be expected to require as inputs a description 
of these gradients. However, because a detailed knowledge of the plasma 
parameters in the aperture region is not generally known before a 
thruster is built and tested, it was considered desirable to require as 
inputs to the model a minimum of plasma parameters (density, potential, 
etc.). This condition was imposed on the model in order that the end 
result be at least reasonably easy to apply in the design phase of a 
thruster. Because it is essentially the difference in these plasma 
properties between the cathode and main discharge regions which drives 
the electrons through the aperture, it was felt that plasma property 
measurements to be used as inputs to the model could be limited to one 
position in the cathode region and one in the main discharge region. 

The plasma density, potential and magnetic field variations 
through the aperture have all been discussed thus far as being impor- 
tant to the processes on going in the aperture region. In addition 
one might expect that the physical area of the aperture as well as the 
magnitude of the current through it might also be important. Thus a 
systematic variation of these quantities suggests itself as a means of 
testing the model . 



III. THEORY 


Development of Theoretical Model 

The approach to the problem of developing a theory which relates 
the current through the baffle aperture to the geometry, magnetic field 
and plasma properties in this region was to treat the electron popula- 
tion as a fluid. Only those forces which act on fluid elements were 
considered; individual particle motions as well as ion motions were 
neglected. The fluid model is appropriate for the following reasons. 
First there are far too many particles involved to follow individual 
particle motions and second the magnetic field can play the role of 

collisions in the sense that it limits electron free streaming and 

22 

forces the electrons to exhibit collective behavior. For motions 

perpendicular to the magnetic field then, the fluid theory is a good 
22 

approximation. The fluid equation of motion for the electrons, known 

23 

as the momentum transfer equation is written, 

m e n [H - + (v • v) v] = -ne(E + v x B) - vP - vnm g v . (12) 

This equation determines the average velocity of the system of particles 
under consideration and assumes that the pressure is locally isotropic. 
Assuming there is no net change in the average fluid velocity and for 
steady state conditions Eq. (12) becomes, 

-V “>* ->- 

ne(E + v x B) + vP + vnm g v = 0 . (13) 

The current density, j, is introduced through the equation, 

j = -nev (14) 

and using the definition of plasma potential in terms of the electric 
field, W = -E, Eq. (13) can be written as, 
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n e(W + - VP + J = 0 (15) 

At this point it is appropriate to specify the coordinate system to be 
used. The coordinate system is chosen such that the axes are everywhere 
parallel or perpendicular to the magnetic field lines. This is il 1 us- 

I 

trated in the cathode pole piece shown in Fig. 5. In this coordinate 

system the current density can be written as, 

-y ~y ~y 

j = Jjl il + j e i 0 + 0.1 in O 6 ) 


and the magnetic field as, 


B = B i„ (17) 


where i a , i e and i lt are unit vectors in the perpendicular, azimuthal 
and parallel directions respectively. Writing Eq. (15) in component 
form yields. 


m v 


Normal to B 

: nev^V - V X P + j Q B + ® 

J'x = 0 

(18) 


m v 



Azimuthal 

: nev Q V - 7„P - J,. B ♦ -f- 

<_J. 

CD 

II 

o 

(19) 


m v 

Parallel to B: nev (J V - v„P + j„ = 0 . (20) 


Because the electron mobility along magnetic field lines is much 
greater than the mobility across the field it is reasonable to assume 
that there are no steady state potential or density gradients parallel 
to the magnetic field. Thus Eq. (20) becomes. 


J II 


= 0 


( 21 ) 


that is, there is no net current flow along magnetic field lines. 
Similarly it can be assumed that no steady state potential or density 
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Fig. 5 Definition of Local Coordinate System 
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gradients exist in the azimuthal direction due to symmetry. Thus 
Eq. (19) can be written as, 


i = — i 
J 9 v 


( 22 ) 


where as = eB/m g is the electron cyclotron frequency. This equation 
says that no current can flow out of the aperture, j x , unless there is 
also a current in the azimuthal direction.^ Finally Eq. (18) can be 
solved for j x to obtain, 

i. - - v ) • 


(23) 


This equation can be written in terms of the classical diffusion co- 

24 

efficient defined by. 


to yield. 


D C • 


j i. = ' D c ( W } (neV ^ V - p ) * 


(24) 


(25) 


For the magnetic field strengths and densities in the aperture region 
it was found that, w/v >> 1, (see Appendix B), therefore v 2 may be 
neglected relative to a> 2 in the denominator of Eq. (24) to obtain, 

(26) 


kTm v 

D c = 


Diffusion coefficients calculated from Eq. (26) have often been found 

12 25 

to be inadequate for use in typical ion thruster plasmas, ’ pre- 
dicting coefficients sometimes orders of magnitude smaller than 


Plasma oscillations leading to potential gradients in the azimuthal 
direction have been proposed as a mechanism for anamolous diffusion. 
However, it is asserted here that no steady state potential gradients 
exist, while maintaining that anamoulous diffusion is still possible. 
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experimentally observed values. This failure of the classical theory 

is attributed to the enhancement of the diffusion mechanism by plasma 
26 

turbulence. A simple and well known semiemperical approach to tur- 

27 

bulent electron diffusion was given by Bohm as. 




_ kT 
" 16eB * 


(27) 


Although this theory has been shown to work much better in the main 

1 2 25 

discharge chamber of an ion thruster than the classical theory, ’ 

27 

it is still only accurate to within a factor of 2 or 3. Calculations 
presented in Appendix B indicate that Bohm diffusion also yields better 
results in the baffle aperture region than classical diffusion theory. 
Therefore in spite of the factor of 2 or 3 uncertainty, Bohm diffusion 
will be used in the remainder of the analysis because of its simplicity 
and relatively good agreement with experimental results. Using Eq. (27) 
in place of the classical diffusion coefficient in Eq. (25) yields. 


Jx = “ i^B (nev ^ V ’ P) 


(28) 


d * 

Now making the change of notation, v x ^ , Eq. (28) can be written 
in the form. 


1 fna dV df\ 
16B (ne dr " dr' 


(29) 


or. 


-l6j x Bdr = endV - dP . 
28 

Integrating over the diffusion depth yields. 


(30) 


This change of notation is appropriate when the net electron current 
flow out of the aperture is in the radial direction. If the current 
flows in the axial direction then v x -* ^ should be used. 
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-16/jj.Bdr = e/ndV - fdP . (31) 

The current density, jj., can be written as, 

= } (32) 

"fc 

where. A, is the area through which the current (I) flows in the 
aperture. Neglecting ionizations in the aperture region the current, 

I, must be constant to satisfy continuity. The area through which the 
current flows however is not constant due to the cylindrical geometry 
and magnetic field configuration, and therefore it must be left under 
the integral sign. Thus Eq. (31) can be written as, 

-161 j| dr = e/ndV - JdP (33) 

or 

I = - e / nd ^ - / dP . (34) 

TS/fdr 

This is the desired theoretical relationship between the current 
through the aperture and the geometry, magnetic field and plasma prop- 
erties around the baffle aperture region. 

Simpl ifications 

Although Eq. (34) provides a theoretical relationship between the 
parameters of interest it is undesirable, in its present form, for the 
following reason. Evaluation of the integrals on the right-hand side 

The current through the aperture due to the motion of the electrons 
is given the symbol (I) in this analysis and is in accordance with 
the standard current convention. Thus the current denoted by this 
symbol is related to the electron current through the aperture 
described in Eq. (5) by, I = -In • 

a 
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of this equation requires a detailed knowledge of the plasma properties 
through the aperture. This is not in keeping with the objectives of 
the investigation specified in Chapter II. Consequently some simplica- 
tion of Eq. (34) is in order. The second integral in the numerator of 
Eq. (34) can be integrated directly to yield, 

/dP = AP (35) 

where, AP, is simply the difference in electron pressure between the 
main and cathode discharge regions. Since the velocity distribution 

i k 

in the cathode region is approximately Maxwellian, the pressure can 
be determined from. 


P = n kT 
c c c 


(36) 


where the subscript "c" refers to cathode region properties. The elec- 
tron pressure on the main discharge side is the sum of the pressures 
due to the Maxwellian and primary electron populations, that is 


P m = n kT + P n 
m mm p 


(37) 


where the subscript "m" refers to main discharge region properties and 

Pp is the pressure due to the primary electrons. This pressure can be 

29 

found from the relation 

P„ = n m w w (38) 

p p e 

where m g is the electron mass and w is the average primary electron 


Some Langmuir probe trace measurements did indicate the presence of 
'v/IO eV primaries in the cathode region. Presumably these primaries 
are created by being accelerated from the cathode orifice across the 
electrostatic sheath which separates the plasma from the cathode, and 
into the cathode region. Inclusion of these primaries in the data 
analysis was found to have a negligible effect, and they will there- 
fore be neglected. 
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speed. This speed can be determined from, 

h m VI 2 = e 


(39) 


where e p is the primary electron energy. Combining Eqs. (35), (36), 
(37), (38) and (39) yields. 


/dP = AP = n kT + 2n e„ - n kT . 
1 mm p p c c 


(40) 


Simplification of the term, JndV, in Eq. (34) is a little less straight- 
forward. Figure 6 shows the density and potential variation through a 

•k 

typical baffle aperture. These and other data not presented suggest 
that a reasonable approximation to the term, /ndV, can be obtained from, 

/ndV = nAV (41) 

where AV is the change in potential across the aperture and n is the 
average value of the density defined by. 


n + n 
n = - m - ~ c 


(42) 


Combining Eqs. (34), (35) and (41), yields. 


I * 


enAV - aP 
1 6/| dr 


(43) 


where AP and n are determined from Eqs. (40) and (42) respectively. 
Equation (43) now is an approximate version of Eq. (34) and requires 
as inputs plasma data from only two points, the cathode and main dis- 
charge regions. The integral in the denominator of Eq. (43) was left 
intact because the strength and configuration of the magnetic field 


These measurements were actually made in the radial field thruster 
described in the next chapter. Tn this thruster the direction 
perpendicular to the magnetic field (horizontal axis) is the radial 
direction. 



ELECTRON DENSITY, n (m 3 ) 
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Fig. 6 Typical Plasma Density and Potential Variations 
Through the Baffle Aperture 


PLASMA POTENTIAL, V (volts) 
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as well as the aperture area are under the direct control of the de- 
signer, and thus the integral can easily be specified, a priori. 



IV. APPARATUS 


Although the aperture region was of prime consideration for this 
study it must be investigated as part of a complete thruster. There- 
fore, the two very different thruster configurations described below 
were tested. 

Radial Field Thruster 

30 

Tests were conducted on the radial field thruster shown schemat- 
ically in Fig. 7. For this thruster and cathode pole piece geometry 
the electron current emitted by the cathode flows axially into the 
cathode discharge region. From there the current is driven radially 
out of the cathode region across the aperture region magnetic field. 

The electrons diffuse across the magnetic field lines in the baffle 
aperture until they reach a field line which carries them into the 
main discharge region. This field line is called the critical magnetic 
field line. 

This thruster consisted of a 14 cm diameter source equipped with 
a set of dished accelerator system grids. The accelerator grid had a 
54 percent open area fraction and was maintained at -500 volts, while 
the screen grid had a 67 percent open area fraction and was operated 
at +1000 volts. Eight radially oriented electromagnets were used to 
control the strength of the magnetic field in the main discharge region. 
The thruster was equipped with the magnetic baffle/pole piece assembly 
shown in detail in Fig. 8. This magnetic baffle assembly allowed the 
magnetic field strength in the aperture region to be varied over a 
wide range of values through the use of a 12 turn solenoid. The 
magnetic coupling rods shown were made of soft iron (as were the pole 
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Fig. 7 Radial Field Thruster Equipped with Magnetic Baffle 
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Fig. 8 Cathode Pole Piece/Baffle Assembly for Radial Field Thruster 
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piece and baffle plate) and were of suffient cross section to assure 
that the magnetic flux through the rods was well below the saturation 
point at the maximum solenoid current used. Since the cathode pole 
piece was part of the main thruster's magnetic circuit, changing the 
current through the radial magnets changed the magnetic field strength 
in the aperture as well as in the main discharge region. However, 
changing the field strength in the aperture by using the baffle mag- 
netic solenoid did not significantly change the magnetic field in the 
main discharge region. This allowed data for testing the model, to be 
collected over a wide range of magnetic field variations. 

Two cylindrical Langmuir probes were used to measure the plasma 
properties on each side of the aperture. The cathode region probe was 
a 0.076 cm diameter wire, 0.076 cm long supported from a quartz tube 
insulator. It was positioned at a radial distance approximately equal 
to the keeper radius at the axial midpoint of the aperture gap. The 
main discharge region probe was a 0.076 cm diameter wire, 0.123 cm 
long again supported from a quartz insulator. The positioning of both 
probes can be seen in Fig. 8. The main discharge region probe could 
also be swept radially through the aperture to facilitate the measure- 
ment of plasma properties as a function of radial position. 

Multi pole Thruster 

31 

Additional tests were conducted on the multipole thruster shown 
schematically in Fig. 9. This was a 15 cm diameter source equipped 
with the same grid set used on the radial field thruster. Thirty-two 
electromagnets wired in series were used to control the strength of 
the magnetic field protecting the anodes. This thruster was outfitted 
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Fig. 9 Multipole Thruster with Large Cathode Region/Maqnetic 
Baffle Assembly 
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with an abnormally large cathode region/magnetic baffle assembly which 
provided a radically different geometry from that of the radial field 
thruster on which to test the model. The direction of electron cur- 
rent flow through the aperture is also indicated in Fig. 9. 

Again two probes similar in construction to the probes used on the 
radial thruster, were used to measure the plasma properties on each 
side of the aperture. Because the magnetic field created by the baffle 
electromagnets was concentrated in the aperture region only, the cathode 
region probe was positioned out at the radius of the aperture and not 
over the keeper as it was for the radial field thruster. The positions 
of both probes can be seen in Fig. 9. The thruster was operated with 

the cathode region/magnetic baffle assembly extended %1 cm into the main 

31 

discharge region as suggested by Longhurst to assure that electrons 
leaving the baffle aperture were not immediately captured by the mag- 
netic fringe field of the inner most anode ring. Finally the cathode 
was electrically isolated from the thruster body so that the net elec- 
tron current emitted by the cathode could be measured directly during 
operation . 

Operation of both thrusters was conducted in the 1 .2 m diameter, 

4.6 m long vacuum facility at the Engineering Research Center of 

★ 

Another cathode region baffle assembly utilizing two hollow cathodes 
was also tested in this multipole thruster. The results of this 
test are given in Appendix C. 

'fc'ic 

The details of this measurement are described in Appendix D. 
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Colorado State University. The facility pressure was maintained below 
1 x 10" 5 torr for all tests. 



V. PROCEDURE 


In order to test the baffle aperture model developed in Chapter III 
experimentally, Eq. (43) was rewritten in the form, 

enAV - AP = - 16l/| dr . (44) 

This equation is now in a form which states that the plasma properties 
on the left-hand side of this equation are a function of the parameters 
on the right-hand side. That is, the plasma properties in Eq. (44) 
should change in accordance with changes in either the current through 
the aperture (I), the aperture region magnetic field strength (B) or 
the area through which the current flows (A). In addition these 
quantities, the current, magnetic field strength and aperture area can 
easily be controlled by the experimenter. 

Consequently the following experimental procedure was implemented 
for the collection of data on the radial field thruster. Before each 
test run the aperture area was set by adjusting the distance between 
the downstream edge of the cathode pole piece and the baffle plate. 

Three different aperture areas were tested. During thruster operation 
the current through the aperture was indirectly varied by changing the 
current through the discharge power supply. In addition the magnetic 
field strengths in the aperture and main discharge regions were varied 
by changing the currents through the radial and baffle electromagnets. 
The magnetic field strength in the main discharge region was varied 
from W to 45 x 10" 4 tesla, while the field in the aperture region 
was varied over the range of 9 to 110 x 10 -4 tesla. These changes re- 
sulted in an order of magnitude variation in the beam current, from 
^ 40 to 485 mA. At each discharge current, magnetic field and aperture 
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area setting tested the following data was collected. Langmuir probe 
traces, from which the plasma properties of interest can be determined, 
were taken in both the cathode and main discharge regions and the dis- 
charge current, discharge voltage, beam current, keeper voltage and 
electromagnet currents were recorded. In addition to the above men- 
tioned variations the cathode and main neutral flow rates were varied 
over the ranges 80-200 and 400-600 mA eq. respectively. The neutral 
flow rate measurements were made approximately every 30 minutes by 
timing the drop a liquid mercury column in a precision bore glass tube. 
After shutdown of the thruster and its subsequent removal from the 
vacuum system the magnetic field strength through the aperture was 
measured with a two-axis gauss meter. The magnetic field strength was 
measured for each setting of radial and baffle electromagnet current at 
which Langmuir probe traces were recorded. 

The procedure used on the radial field thruster was also followed 
for the collection of data on the multipole thruster with the exception 
that the net cathode emission current was also recorded at each oper- 
ating condition tested. A summary of the parameters varied in the 
tests on both thrusters is provided in Table I. 


* 32 

Both Langmuir probes were cleaned by ion bombardment before any 

traces were taken. Probe traces were analyzed using a program de- 

33 

veloped by Beattie. 
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Table I 

Test Conditions 


Test 

No. 

Thruster 

Aperture 

Area 

(10- 4 m 2 ) 

Flow Rate Range 
(mA eq.) 

Aperture 

Current 

[A] 

Aperture Magnetic 
Field Strength 
(10 -4 tesla) 

Cath. 

Main 

1 

Radial 

5.7 

1 53-200 

650 

3. 0-4. 5 

9-111 

2 

Radial 

6.9 

102-158 

475-510 

2. 3-3. 4 

15-107 

3 

Radial 

9.1 

80-100 

406-500 

2. 2-3. 5 

13-100 

4 

Multi pole 

11 .0 

200-250 

400-450 

1 .5-3.2 

10- 33 




VI. EXPERIMENTAL RESULTS AND DISCUSSION 


Baffle Aperture Current 

Measurement of the net cathode emission current during the opera- 
tion of the multipole thruster as well as the discharge current and 
beam current allowed the two methods developed in Chapter II for deter- 
mining the current through the aperture to be compared. Combining 
Eqs. (5) and (11) yields. 



That is, the net cathode emission current should approximately equal 
the discharge current minus the ion flux to the beam and to cathode 
potential surfaces. A comparison of the theory with the experimental 
results is given in Fig. 10. The best fit was obtained when the effec- 
tive screen grid open area ( A s ) was taken to be the entire area covered 
by the grid. This assumes that every ion that encounters the accelera- 
tor system gets extracted into the beam and none recombine on the screen 
grid webbing. Clearly the data of Fig. 10 suggest that the current 

through the baffle aperture can be determined using Eq. (45) and measured 
anode and beam currents. 

Comparison of Model with Experiment 

To test the validity of Eq. (43) experimentally it was found more 
convenient for the presentation of the results to rewrite it in the 
form, 

/f dr = - [e " AV " AP] ‘ (46) 

Thus the value of the integral, dr, computed from Eq. (46) using 
measured plasma data can be compared to the value determined from 
measuring the magnetic field and area variation through the aperture 
directly. A description of the procedure used to measure the value 
of this integral directly is given in Appendix E. 



CALCULATED CURRENT, I D -I lh (— ^ — * (amp) 
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NET CATHODE EMISSION CURRENT, I E (amp) 


Fig. 10 Comparison of Methods for Determining the Baffle 
Aperture Current 
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The results of the comparison of the integrals calculated by the 
theory (Eq. (46)) to the directly measured values for the radial field 
thruster is given in Fig. 11, where the solid line represents a linear 
least squares curve fit forced to go through the origin. This line has 
a correlation coefficient of 0.83 and a slope of y. The different 
symbols represent the different test configurations listed in Table I. 

The theory is seen to yield remarkably consistant results even over ex- 
treme changes in thruster operating conditions. The scatter observed 
in the data is believed to be due at least in part to experimental error. 
The single largest source of error being the analysis of Langmuir probe 
data. The amount of scatter is believed to be good when compared to 
other studies of this nature. Ideally, however, the curve fit should 
have a slope of unity if the model is exactly correct. As it stands, the 
model predicts values of the integral that are a factor of two low. It is 
noted however that this is within the accuracy of the Bohm diffusion co- 
efficient. Consequently a straightforward method to obtain good agree- 
ment between the model and experiment would be to change the constant 
jg in the Bohm diffusion formula. The value yg has no theoretical justi- 
fication but is an empirical number agreeing with most experiments to 

34 

within a factor of two or three. An experimental fit of this constant 

35 26 

by Spitzer yielded a value of 0.21 and Yoshikawa and Rose give the 

constant in the range 0.031 to 0.063 depending on the magnitude of the 

density fluctuation. However, both of these studies as well as Bohm's 

original work were conducted at magnetic field strengths and plasma 

densities much greater than those normally found in an ion thruster. In 

the few diffusion studies actually conducted in operating ion thrus- 

1 2 25 36 36 

ters 5 0 which employed Bohm diffusion only Robinson obtained 



CALCULATED INTEGRAL, fi dr = [en AV-Ap] , [tesla /m] 


i * 


Fig. 
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results which were good to better than a factor of two or three. Con- 
sequently it is not too surprising that the coefficient is off by a 
factor of two. A more important result is the minimum of scatter ob- 
served over substantial changes in operating parameters. Good agreement 
between the model and experiment would be obtained in this case for an 
empirical constant of -g- rather than yg . 

It was of interest to test the model on a radically different 
cathode region and thruster geometry from that of the radial field 
thruster. The multipole thruster of Fig. 9 meets this requirement. 

The major differences between this thruster geometry and that of the 
radial field thruster are the following. First, the current flows es- 
sentially in the axial direction through the aperture rather than 
radially. Second, the large size of the cathode region means that the 
mean diameter of the aperture is considerably larger in this case as is 
the distance between the cathode itself and the aperture. This makes 
the coupling between the cathode and the anode more difficult. Finally 
the main discharge region is essentially magnetic field free making less 
obvious the location of the critical magnetic field line from which 
electrons find access to the main discharge. Due to the axial current 
flow direction through the aperture for this thruster, Eq. (43) must be 
written in the form, 

j|dz = - ygj (enAv - AP) (47) 

g 

where z is in the axial direction. The measurement of J^- dz, for this 
geometry is also discussed in Appendix E. A comparison of results ob- 
tained from the calculations of Eq. (47) and the measured integrals is 
given in Fig. 12. For this case it is seen that the theory predicts 
values that are scattered about the line of perfect correlation with no 



CALCULATED INTEGRAL 
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Fig. 12 Comparison of Calculated Integrals with the Measured 
Values for the Multi pole Thruster 
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change of the empircal constant jg required. In addition the solid 
symbols in Fig. 12 indicate that the theory is independent of whether 
or not the high voltage is on as would be expected. It is not clear, 
however, why the model should agree correctly for this thruster geometry 
and be off by a factor of two for the radial field thruster geometry 
since one would expect the diffusion processes to be similar in both 
cases. It is possible however that a change such as that in the 
cathode pole piece geometry could change the nature of the turbulence 
which affects this empirical constant. It should also be pointed out 
that a considerably smaller range of values for the integral, dz, 
was tested on the multipole thruster. This was primarily a result of 
the poor coupling between the cathode and the anode which prevented 
operation at higher aperture region magnetic field strengths. In any 
case this simple model has been shown to be correct to within a factor 
of 2 for extremely different cathode region geometries. In addition 
for each particular thruster geometry the theory is consistant over 
very substantial changes in operating conditions. 

Inspection of Eq. (43) reveals that the electron current is driven 
through the aperture under the influence of two forces; the force due 
to the potential difference, enAV, and the force due to the pressure 
difference, AP. It was of interest to calculate the relative magni- 
tudes of these forces to determine which has the dominate effect. 

Table II lists the values of these terms for several cases covering a 
wide range of operating conditions. In addition the terms, n c kT c , 

n kT and 2n„e„ which comprise AP according to Eq. (40) are also 
mm p p r ^ 


tabulated. 
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Table II 

Comparison of Driving Forces 
(Percent of Total ) 


Thruster 

Geometry 

enAv 

AP 

n kT 
m 

2n e 
P P 

n kT 

C i 


i%) 

it) 

(*) 

(*) 

it) 

Radial 

90 

10 

- 3 

-2 

15 

Radial 

90 

10 

- 1 

-0.4 

12 

Radial 

94 

6 

- 8 

-6 

20 

Radial 

98 

2 

- 7 

-0.3 

8 

Radial 

106 

-6 

- 8 

-6 

7 

Radial 

103 

-3 

-16 

-1 

15 

Mul tipol e 

99 

1 

-22 

-14 

37 

Mul tipol e 

74 

26 

- 5 

-4 
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Clearly the force due to the potential difference has the dominant 
effect. In fact in some cases the pressure force is seen to act in the 
direction opposite to the flow of electrons. Further, the contribution 
to the pressure from the primary electrons is seen to be generally small 

Design Usefulness 

The design usefulness of the model is limited due to the factor of 
2 or 3 uncertainty in the empirical constant in the Bohm diffusion model 
and the failure of the experimental results to establish firmly a new 
value for this constant of proportionality applicable to any baffle 
aperture geometry. Further, there is presently no model from which the 
cathode region plasma density can be calculated, and an a priori specifi 
cation of this quantity can only be done approximately. The theory is 
still useful, however, in that the value of the required proportionality 
constant appears to be the same for a given thruster geometry over a 
wide range of operating conditions. Even considering the uncertainty 
the model could be used to provide a starting point for the cut and try 
design procedure. A brief description of how to apply the model to aid 
in the design of the baffle aperture region is provided below. 
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In general it would be desirable to be able to specify what open 
area and magnetic field strength are required for the aperture region 
to effect thruster operation at a given condition. This translates 

n 

into the determination of the integral, dr. It is seen that in order 
to calculate this integral from Eq. (46) (rewritten here for convenience), 

Jj *• - - W C-AV - 4P] (48> 

the user must determine, a priori, the value of each of the quantities 

on the right-hand side of this equation. In general the designer of a 

thruster will know from other considerations the values of the beam 

current, discharge current and primary electron energy at which the 

thruster will operate. In addition the main discharge chamber geometry 

will be defined. Once the beam current has been specified the main 

* 

discharge region density, n m , can be found from the Bohm condition for 
a stable sheath which is expressed through the equation 

I. = A n e /kT /m. . (49) 

i b s m nr i ' ' 

The use of Eq. (49) requires the user to supply a reasonable approxi- 
mation for the main discharge region electron temperature, T , which 
can be taken to be ^5 eV. As stated previously the plasma density in 
the cathode region is the most difficult parameter to specify and 
therefore will be a likely source of error. However, typical cathode 
region densities in this study were found to be of the order of 
1 O 1 8 m" 3 . The electron temperature of the cathode region plasma can 
be taken to be approximately T g = 23000°K (2 eV). The density of the 
primary electrons in the main discharge region is also difficult to 
specify. Fortunately the pressure due to the primaries is small, as 
was shown earlier, and perhaps it can be neglected for these 
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calculations. Therefore the pressure difference, AP, can now be 
calculated in accordance with Eq, (40). In addition, the current 
through the aperture can be calculated from Eq. (11). Finally, the 
potential difference across the aperture can be approximated as the 

* 

difference between the discharge voltage, Vp» and the keeper voltage, 

V as ’ 

AV = V D - V k . (50) 

This difference should also be equal to the energy gained by the primary 
electrons after being accelerated through the aperture. Figure 13 shows 
the energy gained by the primaries, which was calculated as Bp - T c 
in units of electron volts, plotted against the difference between the 
discharge and keeper voltages for the radial field thruster. Clearly 
this approximation is appropriate. Thus all the quantities on the 
right-hand side of Eq. (48) can be determined to one degree or another. 

g 

The value of the integral / ^ dr can then be calculated and subsequently 
the baffle aperture region can be constructed to achieve this value. 

Finally it is mentioned in passing that although the model is only 
certain to within a factor of 2 this is perhaps better than the certainty 


* 

This was not true for the multipole thruster's cathode region used 
in this study because the large distance between the cathode and 
baffle aperture was accompanied by an increase in plasma potential 
over this distance. However this should be a reasonable approxima- 
tion for more conventionally sized cathode regions where the distance 
from the cathode to the aperture region is comparable to the distance 
across the aperture region magnetic field. 



PRIMARY ELECTRON ENERGY GAIN, e D -T c (volts 



Fig. 13 Correlation of the Energy Gained by the Primary 

Electrons Through the Aperture with the Difference 
Between the Discharge and Keeper Voltages for the 
Radial Field Thruster 
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to which the plasma parameters required as inputs to the model can be 
specified. Also, it should be pointed out that the theory says nothing 
about where in the thruster's main discharge chamber the primary elec- 
trons should be injected, or where in the cathode discharge region the 
cathode should be located, or even what fraction of the total propellant 
flow rate should go through the cathode, all of which can influence 
thruster performance. 



CONCLUSIONS 


The following conclusions are drawn based on the study of the 
baffle aperture region. 

1) The Bohm diffusion theory predicts values of the diffusion co- 
efficient for the baffle aperture region which agree better with the 
experimentally determined values than does the classical theory. 

2) The classical (col 1 isional ) theory typically predicts values of 
the diffusion coefficient which are at least an order of magnitude to 
small . 

3) The baffle aperture current can be measured directly by measur- 
ing the net cathode emission current. 

4) The magnitude of the current through the baffle aperture can 

be calculated as the difference between the discharge current and the ion 
current which goes to cathode potential surfaces and into the ion beam. 
These calculated values agree well with the measured values. 

5) The simplified theoretical model of the baffle aperture region 
agrees correctly with the experimental results for the multipole 
thruster/cathode discharge region geometry. On the radial field thruster 
the model differs consistantly from the experimental results by a con- 
stant factor of two. 

6) The potential difference across the aperture is the dominant 
force driving the electron current across the baffle aperture magnetic 
field lines into the main discharge region. 

Future work should focus on the development of a model which can 
be used to determine the plasma properties of the cathode discharge 
region. In addition an in depth study of the relationships between the 
electron diffusion mechanism, the plasma turbulence characteristics 
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(frequency, magnitude) and the cathode discharge region geometry might 
be necessary to determine why there appeared to be a factor of two dif- 
ference between the results obtained on the radial and multipole 
thrusters . 
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APPENDIX A 


MAGNITUDE OF NEGLECTED ION CURRENTS 


In the derivation of Eq. (5) the ion current through the baffle 
aperture as well as the ion current to the cathode region walls were 
neglected. To show that this is appropriate the magnitudes of these 
currents will be calculated and compared to the electron current 
through the aperture for the worst cases observed. 

The flux of ions that reach the cathode region walls is given by 
the formula, r 


r i = Vc 


kT 


m. 


( A— 1 ) 

to e 

1 5 

assuming they reach the sheath edge with the Bohm velocity. Assuming 
singly charged ions, the ion correct to the walls then can be found 
from. 


I .= n eA 
i c c 

to 


m_ 


(A-2) 


Equation (A-2) can now be used to calculate the ion current to the 
cathode region walls. The area available for ion recombination in the 
cathode region was about 


A = 2.5 x 10 -4 m . 
c 

For the highest density case the electron density was 

n = 2.1 x 10 18 m -3 
c 

and the electron temperature was 

T = 17 ,400°K (1 .5 eV). 
c 

Substituting these values into Eq. (A-2) yielded an ion current to the 


I . =0.07 amp 

1 w 


walls of 
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which was less than 2.2% of the electron current through the aperture 
(3.2 amp) as calculated by Eq. (11). Equation (A-2) can also be used 
to calculate the ion current through the baffle aperture if it is also 
assumed that the ions reach the aperture with the Bohm velocity. For 
this case the maximum observed ion current was calculated to be 

I.. = 0.055 amp 

cl 

which was 2.54% of the electron current through the aperture. These 
examples represent the worst cases observed. The average ion currents 
were typically closer to 1% the aperture electron current. 



APPENDIX B 


COMPARISON OF CLASSICAL AND BOHM DIFFUSION COEFFICIENTS 


In the course of the development of the baffle aperture model it 
was necessary to determine whether the classical or Bohm diffusion des- 
cribed electron flow through the aperture region more closely. For this 
reason an order of magnitude comparison of the theories with experi- 
mentally determined values of the diffusion coefficient was performed. 

t 

The experimental diffusion coefficients were calculated by taking 


Eq. (25) 



(B-l ) 


and writing it in the form, 
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Jx 
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kT Ar kT Ar 

This equation computes an approximate average value of the diffusion 
coefficient through the aperture. The average density, n, was computed 
according to Eq. (42) and the average temperature, T, was computed from 
the equation. 
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(B-3) 


which roughly takes into account the presence of primary as well as 

Maxwellian electrons in the main discharge chamber. The distance over 

which the change in plasma properties occurs, Ar, was taken to be equal 
* 

to the gap size for all cases and was assumed to be the same for both 


* 


The gap size was taken to be the distance from the downstream end of 
the cathode pole piece to the baffle plate for the radial field 
thruster. 



the potential and density variations. The diffusion coefficients pre- 
dicted by the Bohm diffusion theory were calculated from the equation. 


d = JsL. 

U B 1 6eB 


(B-4) 


where is the magnetic flux density measured at the midpoint of the 
aperture gap. In addition the diffusion coefficients predicted by 
classical diffusion were calculated according to, 


D, 


kTrn v 
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e 2 B 2 
9 


(B-5) 


The collision frequency, v, was calculated considering only electron- 

neutral and electron-ion collisions. The electron-neutral collision 

38 

frequency, v gn , was approximated from, 

= n Q w 
en n n 

where the neutral number density, n n> can be determined from the theory 

39 

developed by Wilbur. The electron-neutral collision cross section 

24 

was approximated from data given by Jahn and the average electron 
velocity w" was given by, 

I O I.T 

(B-6) 
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In addition the electron-ion collision frequency, v , was evaluated 


from the expression given by Holt and Haskel , 


38 
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where A is the non-dimensional plasma shielding distance which can be 
38 

l t 3/2 1 


calculated from,' 
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Table B-I shows the comparison of the diffusion coefficients calculated 
from the above equations for a number of cases covering a wide range 
of conditions. Clearly Bohm diffusion gives better results than the 
classical theory, the latter being typically at least an order of mag- 
nitude too small. In addition Table B-I lists the values of the elec- 
tron cyclotron frequency, w = eB/(2imi e ) and the collision frequency, 
v, at these conditions. In all cases it is seen that the assumption 

oi/v » 1 is satisfied. 

\ 


Table B-I 

Diffusion Coefficients 


Thruster 

B 

g 

O) 

V 

D 

exp 

D B 

D c 

Geometry 

[1 0 -4 tesla] 

[Hz] 

[Hz] 

[m 2 /s] 

[m 2 /s] 

[m 2 /s] 

Radial 

16.7 

'4 . 6x1 0 7 

1.4xl0 6 

200 

TOO 

8.0 

Radial 

92.7 

2.5xl0 8 

8.4x1 0 6 

180 

120 

1.0 

Radial 

58.4 

1 .6x1 0 8 

1.5xl0 6 

43 

37 

0.9 

Radial 

42.1 

1 .2x10® 

2.8x1 0 6 

53 

37 

2.2 

Multi pole 

16.9 

4.8x1 0 7 

1 .7xl0 6 

no 

98 

9.1 

Multi pole 

27.0 

7 .5x1 0 7 

2.4x1 0 6 

26 

66 

5.3 



APPENDIX C 


MULTIPLE HOLLOW CATHODES 

Plasma property measurements were made in a 15 cm multi pole thrus- 
ter equipped with two hollow cathodes in an effort to test the baffle 
aperture model. This thruster configuration is shown schematically in 
Fig. C-l . An assumption used in the derivation of the model was that 
the plasma in the cathode region would be azimuthally symmetric. This 
requirement was easily satisfied for the two cathode region geometries 
described earlier due to symmetry. However for the cylindrical shell 
geometry of this cathode discharge chamber and the azimuthally asymmetric 
location of the cathodes, this requirement is not readily satisfied. 

The first step to attempt to produce the desired symmetric plasma was 
to achieve operation at equal emission currents from both cathodes. 

Toward this end the cathodes were electrically isolated from the thruster 
body, as described in Appendix D, to facilitate measurement of the 
emission currents. Next, some method of controlling the current from 
each cathode was required. This control was effected by placing variable 
resistors in the leads to each cathode in the manner suggested by 
Fig. C-2, where the resistance in each lead was adjusted until the 
emission currents were equal. Control of the emission currents could 

also be achieved by adjustment of the cathode heater powers and neutral 
37 

flow rates. Unfortunately, regardless of the control scheme employed, 
the plasma in the cathode region was not azimuthally uniform even when 
both cathodes were emitting equally. The variation of density with 
azimuthal position (denoted by the angle e from the right cathode) is 
shown for example in Fig. C-3. Clearly the electron density in the 
cathode region (circular symbols) is not uniform, in fact, it even falls 
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Fig. C-l Multipole Thruster Equipped with Two Hollow Cathodes 
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FOR EMISSION CONTROL 


Fig. C-2 Location of Variable Resistors for Control of Relative 
Cathode Emission 



ELECTRON DENSITY (m -3 ) 


t 


XIO ' 7 



Fig. C-3 Aximuthal Plasma Variations 
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below the main discharge region density (square symbols) away from the 
cathodes. In spite of the asymmetry of the plasma in the cathode region 
the main discharge region plasma is seen to be quite uniform. Presum- 
ably, this is a consequence of the low magnetic field strength in this 
region which is characteristic of multipole thrusters. The thruster 
was operated over a wide range of conditions in several unsuccessful 
attempts to achieve a uniform cathode region plasma. Even reducing the 
aperture gap size from 5 mm to 2.5 mm did not facilitate a significant 
improvement in the uniformity. It is postulated that this azimuthal 
asymmetry is a result of the loss of ions to the cathode region walls 
in the azimuthal direction and the subsequent loss of electrons through 
the aperture. That is, the hollow cathodes act essentially as point 
plasma sources and the cathode region walls serve as distributed ion 
sinks. This combination could qualitatively produce the large azimuthal 
density gradients seen in Fig. C-3. 

Finally, regardless of the exact mechanism involved, the proposed 
aperture region model could not be tested on this cathode configuration 
because of the cathode region plasma asymmetry. This asymmetry prohibits 
the one-dimensional ization of the diffusion equation. When plasma prop- 
erty gradients in both the azimuthal and normal directions are included 
in the analysis the equation for the current density through the aperture 
(j x ) determined from Eqs. (18) and (19) becomes. 


v e 


J_L = - 


m g (v' i + a) ) 


2x (nev x V - v x P) + 77-2—2 , (nev 0 V - v 0 p) 


B(v +o) ) 
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and the current density in the azimuthal direction is given by, 
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The presence of the azimuthal plasma property gradients in Eq. (C-l) 
prohibits the use of the simplifications applied in Chapter III to the 
one-dimensional analysis. Consequently, plasma property measurements 
to be used as inputs to these equations could not be limited to one 
point from the main discharge region and one from the cathode region. 
Further, requiring plasma data to be used as inputs to the model from 
more than two locations would make the model prohibitively difficult to 
use as a design tool . 

In addition, it is pointed out that if the Bohm diffusion coeffi- 
cient is incorporated into Eq. (C-l) in the same manner as before and 
neglecting v 2 relative to <d 2 Eq. (C-l) becomes, 

jx = - t^b ( neV * V - + f (neV Q V - V 0 P) . (C-3) 

This equation indicates that even small azimuthal gradients (relative 
to the gradients in the normal direction) can have a significant effect 
on the current through the aperture. 



APPENDIX D 


ELECTRICAL ISOLATION OF CATHODE 

In order to measure the emission current from a hollow cathode in 
an operating ion thruster it is necessary to isolate electrically from 
the thruster body. In this study the cathode was isolated from the body 
with boron nitride insulators and electrical connections were made in 
accordance with Fig. D-l. These connections included the addition of 
a lead from the negative side of the anode power supply to the thruster 
body to supply electrons to neutralize the ions reaching the body. The 
location of the ammeter shown in this figure allows it to measure only 
the net cathode emission current, I^, (i.e., the current emitted by the 
cathode and not subsequently collected by the keeper). The extension 
to additional cathodes is straightforward, where additional keeper, 
vaporizer, and cathode heater power supplies are required for each 


additional cathode. 




Electron Current to Thruster Body 


Fig. D-l Electrically Isolated Cathode Circuit Diagram 




APPENDIX E 


APERTURE REGION MAGNETIC FIELD MEASUREMENTS 

In order to test the validity of Eq. (46) it was necessary to de- 
termine the values of the integral on the left-hand side of this 
equation by directly measuring the area through which the current flows 
and the magnetic field variation through the aperture. A closeup 
drawing of the radial field's cathode pole piece baffle aperture region 
is given in Fig. E-la. For this particular geometry it is assumed the 
electron current flows through increasingly larger diameter right 
circular cylinders as it traverses the baffle aperture region as sug- 
gested by the dashed lines in Fig. E-la. In reality the current flow 
is normal to the magnetic field vectors 17 but the error introduced by 
using the area variation of Fig. E-la is believed to be negligible 
since the magnetic field vectors through the aperture point primarily 
in the axial direction. Figure E-lb shows the normalized area varia- 
tion through the aperture as determined from Fig. E-la. The area was 
normalized according to the formula, 

A* = A/A g . (E-l ) 

In addition it was found convenient for data analysis purposes to 
normalize the magnetic field strength and radial position according' 
to the formulas, 

B* = B/G g (E-2) 

and 

r* = r/r g (E-3) 

where the stars indicate dimensionless quantities and the subscript "g" 
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Fig. E-la Radial Field Cathode Pole Piece Baffle 


Aperture Region Showing Area Through which 
the Electron Current is Assumed to Flow 



NORMALIZED RADIAL POSITION, r* 

Fig. E-lb Example of Normalized Area Variation Through the 
Aperture Determined According to Fig. E-la. 
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refers to quantities pertaining to the center of the aperture gap. 

The magnetic field variation through the aperture was measured at ten 
radial loations with a two axis gauss meter. This was done for each 
gap size and current setting of the electromagnets tested because the 
shape as well as the strength of the magnetic field in this region de- 
pended on the values of these parameters. A typical example of the 
magnetic field variation normalized with the magnetic field strength 
in the gap (B^) is given as a function of the non-dimensional ized radial 

position (r/r ) in Fig* E-2. If at each radial position the value of 

* * 

the magnetic field strength (B ) is divided by the area (A ) from 

★ •k 

Figs. E-2 and E-lb respectively, the result is the variation of (B /A ) 
through the aperture given in Fig. E-3. Now integration of this curve 

k k k 

will yield the value of the non-dimensional ized integral /B /A dr . 

The distance over which this integration is performed is of primary 
importance. The proper value of the integral is determined by integ- 
rating over the distance between the radial location of the cathode 
region probe and the critical magnetic field line. Using Eqs. (E-l), 
(E-2) and (E-3) the dimensional value of the integral can be recovered 
from, 

n B r * 

/ f d r = d r . (E-4 

"g A 

As mentioned previously, this was done for each different electro- 
magnet current setting and aperture area. The same procedure was 
followed for the determination of the integral, / ^ dz, for the multi- 
pole thruster configuration. It was noted in this case, however, that 
the field shape through the aperture was independent of the current 
through the main electromagnets. This is because the baffle aperture 
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0.0 0.2 0.4 0.6 0.8 10 1.2 


NORMALIZED RADIAL POSITION, r*,r/r g 

Fig. E-2 Example of Normalized Magnetic Field Strength 
Variation Through the Aperture 
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Fig. E-3 Example of Variation of the Ratio B*/A* 
Through the Aperture as Determined from 
Fig. E-l b and E-2 ' 
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pole piece was not part of the main discharge region's magnetic circuit. 
In addition, evaluation of the integral was performed over the dis- 
tance between the two probes due to the vagueness of the location of 
the critical magnetic field line in this thruster. 
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